Topological Field Theory for p-wave Superconductors 
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We propose a topological field theory for a spin-less two-dimensional chiral superconductor that 
contains fundamental Majorana fields. Due to a fermionic gauge symmetry, the Majorana modes 
survive as dynamical degrees of freedom only at magnetic vortex cores, and on edges. We argue 
that these modes have the topological properties pertinent to a p-wave superconductor including the 
non-abelian braiding statistics, and support this claim by calculating the ground state degeneracy 
on a torus. We also briefly discuss the connection to the Moore-Read Pfaffian quantum Hall state, 
and extensions to the spinful case and to three-dimensonal topological superconductors. 

PACS numbers: 74.20.De, 74.20.Rp, 03.65.Vf 



Introduction - The great current interest in topological 
phases of matter is mainly due to the intrinsic scientific 
challenges to discover, engineer and describe them. The 
hallmark of a topological phase is that the elementary ex- 
citations carry unusual quantum numbers, and that the 
quantum mechanical wave functions acquire phase fac- 
tors that are insensitive to noise, when these excitations 
are braided around each other. It has been proposed to 
code quantum information in these phase factors, with 
the prospect to greatly improve both quantum memo- 
ries and quantum processors this possible technologi- 
cal application is fascinating. The field was initiated by 
the discovery of the quantum Hall (QH) effect, fueled by 
the theoretical proposal for topological spin liquids, and 
recently by the surge of investigations of topological in- 
sulators and topological superconductors (SC)[2j]. The 
theoretical study of these systems range from detailed 
microscopic models based on realistic band structure, to 
general descriptions in terms of effective low energy theo- 
ries. An important class of the latter are the topological 
field theories (TFTs) that distill the topological informa- 
tion about the states of matter. Prime example of such 
TFTs are the Chern-Simons (CS) theories that describe 
the incompressible quantum Hall fluids. 

The perhaps most familiar system that is a topolog- 
ical phase is an ordinary s-wave SC, where the non- 
trivial braiding phase is the minus sign picked up by the 
wave function as a spinon (a neutral spin half Bogoli- 
ubov quasiparticle) encircles a vortex of strength h/2e 
in the Cooper pair condensate Q. That this is a bona 
fide topological phase was stressed by WenQ, and it was 
later realized that the corresponding topological fleld the- 
ory is a BF theory, which in two space dimensions is 
closely related to the CS theories For the case of 

three space dimensions it was known that the BF theories 
could provide a topological mass 7], and the connection 
to superconductivity was discussed in severalpapersQ. 
A BF theory for d-wave SC was conjecturedl6|, and de- 
rived and analysed by Hermanns [gi]. Recently it has also 
been shown that the TFTs describing topological insula- 
tors are of the BF typeflo']. It is thus natural to assume 



that also a chiral p-wave superconductor should be de- 
scribed by a BF theory, but it is also obvious that such 
a theory must be more complicated in order to correctly 
describe the gapless Majorana modes known to exist both 
in vortex cores, and at boundaries [11., .12J. Furthermore, 
the vortices must satisfy the braiding rules appropriate 
for Ising type non-abelian statistics |12h15| . In this letter 
we propose a TFT for a two-dimensional spin-less chi- 
ral SC that satisfies these requirements, comment upon 
the relation to the Moore-Read Pfaffian QH state, and 
briefly touch upon how to incorporate spin and how to 
generalize to the three-dimensional case. 

Our theory differ from previously proposed effective 
theories for the Moore-Read state in that it has a fun- 
damental Majorana fermion fleld that enters at the same 
level as the topological gauge fields, and just as these it 
does not have any bulk dynamics. In the same way as 
a Chern-Simons gauge field, it does have dynamics on 
edges, but also at vortex cores, which for many purposes 
can be modeled as internal boundaries. The rationale 
for this picture is that it explicitly exhibits the full low 
energy content of the theory. This is in contrast to the 
approach based on non-abelian Chern-Simons theory fl6|, 
where the presence of Majorana zero modes enters implic- 
itly by the choice of representation of the Wilson loops 
that describe the worldlines of the fundamental quasipar- 
ticles. 

The model - Our model is defined by 
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where the second line is obtained by partial integration, 
jv and jq are vortex and quasiparticle currents respec- 
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tively, a and b are gauge potentials, and 7 is a (one- 
component) Majorana fermion. £gf is a gauge fixing 
term that is discussed below. Including only the first 
term in the parenthesis, we retain the BF action, ^ bda, 
for an s-wave SC. Under parity, P, (x,?/) — ^ (— a;, y), 
the fields in ([2]) transform as (ao, ax, ay) -> (ao, — a^, ay), 
{bo,bx,by) -)• (-6o,''a:, and 7^7 and un- 
der time reversal, T, as (ao,aa;,aj^) — > (ag, — a^:, — Oj,), 
(&o, ^x, &y) — > (—^0, bx,by) and 7^7. As is well known, 
this implies that 6(ia is invariant under both P and T, 
while the new term ijd'yda violates both these symme- 
tries as appropriate for a chiral SC. 

We now give a heuristic argument, which will be made 
more precise later, for why ^ and ^ describe a SC 
with dynamical Majorana modes localized on vortices 
and edges. Up to boundary terms, ^ is invariant un- 
der the two gauge transformations, — )■ -(- S^jAq and 
bfi ^ b^ + d^Ab, and it is known that in the s-wave case 
the boundary supports bosonic degrees of freedom that 
will generically be gappedQ. In the present case, this 
is not true as will be discussed later. The terms in the 
second line of ^ possess another local symmetry, 

7 ^ 1 + 29, (3) 

^ &p - - tidf,e/2 - iOdf^e , 

where is a Majorana fermion field. This symmetry is 
violated by the terms on the third line of ©, i.e. by 
the total derivative term and the one proportional to the 
vortex current. This implies that the Majorana fermion 
can be gauged away everywhere except at the boundary 
and at vortex cores, where it emerges as a dynamical 
degree of freedom, in close an alog y to the bosonic edge 
modes in CS or BF theories [gI Il7l| Using the suggestive 
notation 

&^ = z75^7/4, (4) 

we note that ([2]) is quadratic in a, b and b, and linear in 
derivatives. 

Integrating the Lagrange multiplier field b, yields da — 
TTjv, showing that a unit strength vortex carries a flux tt. 
Using the notation ^ for the pertinent Green's function, 
we have a = ^jv + dAa and substituting this in ^ gives 
the bulk Lagrangian 

1 TT 

^ = ^Jv 7*9^7 - jqAp = -^maj - jv Jq , (5) 

where we ignored the gauge-fixing term. The last term 
in (O is an abelian statistical interaction that gives a tt 
phase for a vortex circling a quasiparticle. If there were 
no internal structure to the vortices, i.e. no Majorana 
modes, two vortices with the same charge would be iden- 
tical bosons. 

To exhibit the Majorana modes, we consider a vor- 
tex source corresponding to N Wilson loops Wc^ = 
exp(z/^^ dxf^b^), 

Ji^i^n = E / drSH^^ - <{r))^ , (6) 

a=l -^0 



where the loop Cq is parametrized as x^{t). Substituting 
in ([5]) yields the action 

'5maj = TE/ dr -faiT)idrJaiT) , (7) 
a=l 

where 7a (t) = j{x^{t)). For the simple case of the loops 
extending over all times, and crossing any fixed time sur- 
face only twice, we can use the time, t, as a parameter 
to rewrite this as, 

'5maj = TE/ dt-1a{t)idaa{t), (8) 

which, as promised, describes a collection of gapless Ma- 
jorana fermions living on the world lines of the vortices. 

Braiding and statistics - To demonstrate that ([S]) re- 
produces the Ising braiding statistics of the vortices, we 
switch to an operator description, where (|S]) implies the 
canonical commutation relations, 

{%{t),%{t)}^28ab. (9) 

Since there is an even number of vortices we can form the 
complex Dirac fields iIja = (72yi-i + «72a)/2, where A = 
1,2...A^, which together with -0^ satisfy = 

5ab and {?/'a,V'b} = {V'a'V'L} = 0- The pairing of 
Majorana fermions in the Dirac fields is arbitrary, and 
amounts to a choice of basis in the Hilbert space spanned 
by the vectors \ai---aN) = \oii) ®---® \oiN), where 
= 0, 1 and = |1) and V'aII) = |0). The 

full description of the low-energy sector, in the absence 
of quasiparticles, is thus in terms of the state vectors, 
\xi . . .X2N',oti ■ ■ ■ apf), since the sources determine the 
gauge potentials a{x,t) and b{x,t) up to gauge transfor- 
mations. Also note that the Hamiltonian corresponding 
to the action ([5]) is identically zero, so the time evolu- 
tion is entirely determined by Berry matrices. From gen- 
eral principles, we also know that these statistical ma- 
trices will only depend on the braiding patterns, since 
they will constitute a (non-abelian) representation of the 
braid group. This group is generated by the elementary 
exchanges Tq : Xa Xa+i, that do not involve any braid- 
ing with the other vortices, and these exchange processes 
can be viewed as time evolutions, 

U{T)\xi ... fa, Xa+l . ■ . X2N;ai " ' • aAr)o(lO) 

= E ^<---"'n/2^^^ ■ ■■^a+l,Xa . ..X2N;a[ ■ ■ ■ a'pf)T ■ 
{a'} 

where the unitary operator U (T) interchanges the posi- 
tion of the a-th and (a -I- l)-st vortex, leaving the oth- 
ers unchanged. We assume that the Majorana oper- 
ators, 7a, which are attached to the vortices, are in- 
terchanged when the corresponding vortices are inter- 
changed, U{T)%UHT) = Ca7a+i and U{T)%+iW{T) = 
C,a+ila- Since 7J = 7a, and the commutations relations 
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should be preserved, we must have = il, but any com- 
bination of these signs is allowed. There are, however, 
only two unitary bosonicflEl operators that exchange 7a 
and 7a+i, while leaving all other 7s unchanged. These 
are U± = e^TT'»T'"+i ^ corresponding to Ca = — Ca+i — Tl, 
which are precisely the transformations that imply Ising 
non-abelian statistics as shown by Wilczek and Nayak 
in the context of the Moore- Read QH state[l3|, and by 
Ivanov in the context of a Bogoliubov-de Gennes descrip- 
tion of a p-wave SC[15|. If V'a = (7a + *7a-i-i)/2, then 
U± = e±*f implying e.g. U+\0) = e^f |0) and 

U+\l} = e~*^|l}, where ^\'\lia\n) — n\n). That U is di- 
agonal in this basis follows already from conservation of 
fermionic parity, so the non-trivial result giving the non- 
abelian statistics is the values of the phases. In Ivanov's 
treatment, the crucial relative sign, C,a = — Ca+i orig- 
inates from the coupling of the Bogoliubov-de Gennes 
quasiparticle to the superconducting order parameter, 
and we believe that a similar interpretation is valid also 
here, since the term ^ adb in the topological action iden- 
tifies e^^dibj as a Z2 quasiparticle charge. The above 
argument for non-abelian statistics is based solely on the 
general properties of U± , and it would certainly be inter- 
esting to derive the expression for U± directly by consid- 
ering time evolution in the topological field theory. By 
quantizing the full field theory using canonical methods, 
one could determine the time evolution of the Majorana 
fields 7a, and from there the time evolution of the state 
vectors. This is, however, quite complicated due to the 
singular nature of the constraints. 

Ground state degeneracy on a torus - To further 
strengthen the case for non-abelian statistics we now 
show that, using a natural lattice regularization, there 
are three degenerate bosonic ground states on the torus. 

Following the strategy of Refs. [1] and [l^, we 
solve the constraints in the absence of sources as, — 
diKa + Oil Li, Pi ^ bi + bi = diAp + Pi/L^, where Aa/p 
are periodic and a^, /3j are constant in space; Li is the 
length of the torus in direction i, so = dx^ai and 

Pi — dx^Pi, where the contour Ci winds the torus 
once in the i-th direction, but is otherwise arbitrary. 
Note that Pi is invariant under the fermionic gauge trans- 
formation. We get the Lagrangian L as ttL = e^^aiPj, 

which implies the commutation relations, [a^, Pj] — ine^^ . 
Next we define the operators, 



jAi — 6 



(11) 



^Ibi + i fc- dx^-ydij] 



leCi J-yec, 



with bi — dx''bi, which satisfy {Ax, By} — {Ay, Bx} = 
0, since and Pj obey the same commutation relations 
as ai and bj in the s-wave case. We would thus expect 
that, just as in that case, there would be a 4-fold degen- 
eracy of the ground state as a direct consequence of the 
commutator algebra. The situation here is however more 
subtle, since Pi depends on the fermionic variable and 



the fermionic integration in Bi is necessary to obtain an 
operator that acts only in the bosonic part of the Hilbert 
space only. A direct calculation shows that for Bi to be 
non-zero, we must use anti-periodic boundary conditions 
on 7 in the i-th direction, which is consistent with 7 be- 
ing fermionic. That we need the fermionic integration in 
can also be understood from the Lagrangian point 
of view. If we try to calculate a correlation function in- 
cluding the naive operator e'^' , we notice that this is not 
gauge invariant, and we cannot gaugefix the Majorana 
field on the contour Ci. We can thus naturally associate 
this integral with the operator, and to have a gauge in- 
variant object we must also add the bi contribution to 
get Bi. 

The new complication occurs when we attempt to ex- 
plicitly construct the four ground states. Since Bi mea- 
sures the Oj flux, one would naively expect that the 
four degenerate ground states, just as in the s-wave 
case, should be labeled by thse fluxes. Denoting the Ai 
eigenstates by \sx,Sy) we would have | + ,+), | + , — ) = 
By\+,+), |-,+) = Bx\+,+) and |-,-) = BxBy\ + ,+). 
In our case, this is no longer true. The state | — ,— ) 
is absent since the operator BxBy vanishes due to the 
fermionic integration, and we are left with three ground 
states as expected if the excitations obey non-abelian 
statistics. Heuristically, this "blocking mechanism" 
comes about since the two contours Cx and Cy necessar- 
ily cross, which amounts to having two fermions at the 
same point. A more rigorous argument is given below. 

Discrete formulation and fermionic gauge fixing - To 
prove the above statements about the operators Bi, and 
also to explain the fermionic gauge-fixing, we extend 
the discretized BF theory due to Adams [19| to our La- 
grangian (21). In Adams' construction, the gauge poten- 
tials are real-valued one-forms, a and b, defined on the 
links of two dual cubic lattices K and K. The action is 



1 



d^xb-k^ d^a 



(12) 



where -k^ is a duality operator that relates a p-form on K 
to a (3 — p)-form on K, i.e. a link on if to a plaquette on 
K etc.. The coboundary operator, d^ , is the discretized 
version of the outer derivative d, so one can thus think 
of as a sum of products of a link on the K lattice 
with the dual plaquette it pierces. Adams showed that 
p2l) is invariant under a ^ a + d^ Aa and 5 — > 6 + d^ At, 



and that it reproduces the expectation values of Wilson 
loops as calculated in the continuum theory. Using stan- 
dard arguments, this also implies that the commutation 
relations among the operators Ai and Bj are the same as 
in the continuum. Also, the coboundary operator, d^ , 
satisfies the Liebniz rule, so partial integrations can be 
performed as usual [20|. 

The vortex Wilson loop Wc has the discrete counter- 
part exp(i ^^ -^^gp 5ij), where C is a loop on the dual 

lattice K, and 6^ the gauge one-form on the link (ij). 
We can extend this construction to give a discretized 
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version of our action by simply making the replacement 
where 7^ is a Majorana fermionic 
0-form defined on the sites of the dual lattice K. The 
discrete version of the fermionic gauge symmetry easily 
follows. 

It is now clear that we have to gauge fix all the ji, 
except those living on the loops, and the simplest pro- 
cedure is to introduce the factor Gi'Ji) = Yii^uc '^(^»)' 
where S{^i) = ji is a fermionic delta function, in the 
path integral measure. 

With this we have given a precise prescription for how 
to introduce point-like vortices, and shown that they will 
support dynamical Majorana degrees of freedom. 

The discrete form of the operators Bi directly follows, 

where we note that the field in the exponent is the gauge 
invariant combination j3ij defined above. To evaluate the 
fermionic integrals, we first assume that there is an even 
number of sites in both the x and y directions. It is 
then clear that both B'^ and By are non-zero, since all 
sites on the loop can be covered once, and only once, by 
expanding the exponential. The proper discretized form 
of the product operator is. 

This operator vanishes since there is an odd number of 
fermionic integrals and expanding the exponential will 
only bring down an even number of 7s. Similarly, for an 
even - odd lattice, B'^ and {BxBy)' will give new states, 
while By vanishes, etc.. It is also easy to see that if we al- 
low for fermionic states, corresponding to an odd number 
of vortices, there will be an extra ground state 

Edge states - We already stressed that the total deriva- 
tive term in ^ violates the fermionic gauge symmetry 
and we can directly read off the action for the dynamical 
Majorana edge mode, 

5'od = ^ J dxdt-f{x,t)[a.j;do - aoda;]j{x,t) , (15) 

where x is the coordinate along the edge. The line in- 
tegral f dx ax measures the total vorticity in the system 
which can have both a bulk and an edge contribution. In 
a p-wave superconductor an edge vorticity, corresponding 
to an edge current, is indeed expected. In our topologi- 
cal description this amounts to taking ax as an unknown 
constant. The value of aq is formally a gauge choice, but 
this is somewhat misleading since it in fact determines 
the velocity v = a^jax of the edge modes and thus is re- 
lated to the edge potential. The situation is very similar 
to that for the edge modes of QH liquids, as discussed 
by Wen[l3|. For a circular droplet with radius i?, the en- 
ergy of the fermionic edge modes is E = vl/R where I is 
the angular momentum, which is integer or half-integer 



valued depending on whether the boundary conditions 
are periodic or anti-periodic. For consistency, we must 
thus take periodic boundary conditions if there is an odd 
number of vortices in the bulk, and antiperiodic if the 
number is even. In this way the zero angular momentum 
edge mode will provide the missing Majorana zero mode, 
which combined with the localized bulk Majorana modes 
in a vortex form Dirac fermions. In a careful treatment 
of the gauge fixing, these boundary conditions are likely 
to be necessary for the fermionic path integral not to 
vanish. 

External fields and the Moore-Read QH state - Cou- 
pling to an external electromagnetic potential. A, is ob- 
tained by shifting a — a + ^ as explained in Ref. 0- 
From the work of Reed and Green, we expect that our 
theory should also describe the non-abelian Moore- Read, 
or Pfaffian, QH state which can be viewed as a p-wave 
paired state of composite fermions [ill [2]| . 

A crucial difference between this state and the p-wave 
SC is that flux and charge cannot be separated, so the 
charge e/4 qusiparticle current couples to both a and 6, 

-Cmr = -e^^^l^ + fepRK + ^Ap) (16) 

TT Z 

+ ^e^-Pa^d.A.-fi^a. + ^b^), 

where we also added the coupling to A and omitted the 
total derivative and the gauge-fixing terms. To get a 
physical motivations for the coupling terms in (fT6| . we 
note that the bosonic part can be written on the standard 
form [it}, 

l-e^^-Pticld^Ap ~ lj,rci 

where ci = a, = 6, and t = (1/2,1),!= (1/4,1/2). We 
see that the Ij vector is chosen to describe a Moore-Read 
quasiparticle which is a half vortex with charge e/4. The 
second component in the charge vector tj just expresses 
that Af^ couples to the charge current of the condensate, 
^ db, with unit strength, while the factor 1/2 in the first 
component expresses the flux charge connection in the 
quantum Hall fluid. The unit of vorticity in the vortex 
condensate, ~ db, is a unit (or Laughlin) vortex whcih 
carries a charge given by the filling fraction = 1/2. 
Following the same steps as took us from ((!]) to ([5]), 

-Cmr = Anaj - ^^dA + ^jA ~ ij^j , (17) 

where £inaj differs from that defined in ([S]) by a factor 
of 2 and the presence of a coupling ^ bdA. The second 
term in (|17|) gives the correct Hall conductivity, the third 
term shows that the quasiparticles have charge e/4, and 
the last term provides them with the abelian statistical 
angle 7r/8, which in the conformal field theory descrip- 
tion originates from the charge sector Thus the La- 
grangian (|16p encodes all the pertinent properties of the 
MR state. 
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We can add local source term with the 1 vector Ilk = 
(1/2,1) to describe the charge e/2, semionic, Laughlin 
holes, but it is less clear how to incorporate the neutral 
fermions except as composites of the non-abelian quasi- 
particles. Finally, we note that the coupling of the cur- 
rent j to both a and b has a natural counterpart in the 
discretized theory [TP]. 

Generalizations ~ The generalization of ([l} and ^ 
to include spin, and to T-invariant topological SCs such 
as the proposed 2D p-wave Rashba non-centrosymmetric 
SCfl^], can be obtained by combining several copies of 
as will be described in future workj23|- A 3D s- 
wave SC is described by a TFT of point-like quasipar- 
ticles and string-like vortices [sj, and a natural general- 
ization for a 3D T-invariant topological SC is, S'top = 
J d^x [Cbf - j^pft^ ^ J^bf,,^] , where 

Cbf = -£^""^(6^. + h'^)d.a^ ■ (18) 

TT 

Here the antisymmetric tensor gauge potential 6 couples 
to the world sheets of the strings defined by the source 
•^v^SB- The new term involves the potential 6^^, which 
is defined in analogy with 1^ but this time in terms of 
two Majorana fields 71 and 72 as, 

b,,^ = ^7'^CT^9^7 , (19) 

where 7"^ = (71,72), and cr^ = (1,(7^), with cr^ the Pauli 
matrices. Analogously to the 2D case, there is a local 
fermionic symmetry which implies that the Majorana 
fermions live only on the world sheets of the strings, and 
for the simplest case of a stationary, circular string with 
radius R parametrized by the angle ip, we get the action, 

S-string ^^JdtJ^ d^j^ (^la^dt - 7 (20) 

with CTip = cos ipay — sin ipax ■ This action describes two 
ID Majorana fields living on the vortex. 

Summary and outlook - We have proposed a topo- 
logical action, involving Majorana fields, for a two- 
dimensional spin-less chiral SC that contains funda- 
mental Majorana fields. This action possesses a local 
fermionic symmetry such that the Majorana fields are 
pure gauge degrees of freedom except at vortex cores 
and on boundaries. We have identified these dynamical 
Majorana degrees of freedom and shown that they repro- 
duce the known topological features of chiral SCs, includ- 
ing the non-abelian braiding statistics, and the three-fold 



ground state degeneracy on the torus. We also gave a dis- 
crete version of our model, where the concept of vortex 
core is sharply defined, and where the necessary gauge 
fixing procedure could be unambiguously defined, thus 
justifying the more formal treatment in the continuum 
model. 

There are several interesting directions to pursue. We 
already mentioned the problem of canonically quantizing 
the theory, where one should be able to explicitly cal- 
culate the non-abelian braiding phases, the extension to 
the spin-ful case where the details should be worked out, 
and also the problem of how to consistently include an 
edge vorticity to describe the expected edge currents. 

In addition one should clarify the relation between 
our effective theory by for the Moore-Read state 
QH state, and the low-energy theory based on a SU (2)2 
Chern-Simons thory previoulsy proposed in Ref. jlGj. 
Also, using the discrete model, it should be possible to 
calculate the expectation values of linked vortex loops, 
from which one can extract the non-abelian statistics as 
is done in fisj. A natural way to make the continuum 
theory ^ less singular is to add a (non-topological) 
Maxwell terms for the gauge fields. As explained in 
Q, this gives a finite London length and a dynamic 
plasmon mode. We expect that this will introduce 
gapped fermion modes in the vortices, but that a single 
gapless Majorana mode remains; this remains to be 
shown. 

Although the fermionic gauge symmetry ^ is not a 
supersymmetry, it is nevertheless interesting to specu- 
late that there might be a connection to supersymmetric 
models. An obvious candidate to study would be the su- 
persymmetric extension of the relativistic abelian Higgs 
model that gives a microscopic derivation of the BF the- 
ory for s-wave SCs [1, Q . 

Another challenging task is to investigate the 3D mod- 
els of the type ((T8)) . especially concerning edge modes. 
Also, since the Majorana modes live on strings, there 
might be a connection between the lattice version of 
and the Kitaev chain(2^], which is known to support Ma- 
jorana modes at the ends of open strings. 
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